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Abstract 

The free Bose gas with attractive boundary conditions is an interesting toy model 
for the study of Bose-Einstein Condensation (BEC), because one has BEC already 
in one dimension. Here we study for the first time the imperfect Bose gas with these 
boundary conditions and prove rigorously the occurence of condensation. 
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1 Introduction 

For the free Bose gas with Dirichlet and Neumann boundary conditions, Bose-Einstein Con¬ 
densation (BEC) is rigorously treated in [4,5]. For the mean-held Bose gas with periodic 
boundary conditions BEC is rigorously proved and a detailed analysis of the thermody¬ 
namic limit is given in [3]. The proof is based on bounds on the correlation functions for 
equilibrium states, given in terms of the correlation inequalities [1,2]. The subtle point in 
this proof is the analysis of the singularity around the zero-mode. 

If one considers attractive boundary conditions instead of periodic boundary conditions, 
the problem changes drastically. The free Bose gas with attractive boundary conditions 
is extensively studied in [4,6]. Due to the gap in the one-dimensional one-body problem, 
one has Bose-Einstein Condensation in all dimensions v > 1. An important result is the 
fact that the condensation is a surface effect. In [7] it is computed that the condensate is 
localized at a logarithmic distance from the boundary. 

The subject of this note is to proceed with the imperfect Bose gas with attractive boun¬ 
dary conditions, i.e. the free Bose gas with attractive boundary conditions plus a mean 


^Email: lieselot. vandevenneOf ys . kuleuven .ac.be 
^Email: andre . verbeureSf ys . kuleuven .ac.be 



field term. The first problem that occurs is to express this mean held term in momentum 
space diagonalizing the kinetic energy (the free Bose gas part). As the spectrum of the 
latter one has two strictly negative eigenvalues, say 6^(0) and 6^(1), separated from the 
rest of the spectrum, then one can discuss the corresponding number operators No and Ni 
as being added to the total number operator in the interaction or not. In section 2.3 we 
argue why they should not be present. The argument is essentially based on the fact that 
we want a space homogeneous mean held term. The model is dehned in 2.2. 

In section 3 we give a completely rigorous proof of the occurence of Bose-Einstein Con¬ 
densation for the imperfect Bose gas with attractive boundary conditions. We perform all 
details only in dimension u = 1. From dimension u >2 on, the proof becomes technically 
more tedious. In particular, because of the fact that the condensate is located near the 
boundaries, for higher dimensions the thermodynamic limit for hypercubic boxes is not 
very suitable nor realistic and should in stead be taken with increasing absorbing balls. 
But we leave this extra exercise for a later occasion in which we consider the problem of 
the shape-dependence. 

In the one-dimensional case we remark that the condensation is equally distributed over 
the two negative energy levels. The condensation is localized in the same area as for the 
free Bose gas with attractive boundary conditions, see [6]. 


2 The Model 


2.1 Attractive Boundary conditions 


If one considers a free gas of bosons in an interval [—L/2, L/2] of length L, then the energy 
levels are determined by the one-dimensional Schrodinger equation (with units ^ = 1) 

-A0 = ez,0^. 


with boundary conditions: 
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where a < 0. 

If one considers these attractive boundary conditions, the spectrum consists of two negative 
eigenvalues tending to the same limit —(when L —> cx)) and an inhnite number of positive 
eigenvalues (for L\(j\ > 2): eiik) for k = 0,1,2 ,where 


6^(0) < 6^(1) < 0 < ei,(2) < 61 ,( 3 ) < ..., 

eL(0) = -a2-O(e-^H), 
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The corresponding eigenfunctions are given by 
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for k even, 

for k odd. 


2.2 Hamiltonian 

We consider a one-dimensional system of identical bosons on an interval [—f, f] C IR with 
attractive boundary conditions. The model is specihed by the local Hamiltonians H^j^p 
on the boson Fock space 


Tja _ rp<j 

^L,MF — -^l 
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2 L 


( 2 . 2 ) 


where is the kinetic energy operator with the eiik) the eigenvalues (2.1) of the free 
Laplacian with attractive boundary conditions: 


T[ = '^eL{k)alak 

fcSN 


The operators al = a*{(l)\) and are the Bose creation and annihilation operators 

with the testfunctions the above eigenfunctions of the free Laplacian with attractive 
boundary conditions. The total particle number operator is denoted by Nl = X)a:gn = 
particle number operator corresponding to the positive spectrum by 
Nl- 

OO 

Nl = ^ alttk 
k=2 

We consider a positive coupling constant A G IR'*' for the sake of thermodynamic stability, 
see [3]. 
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2.3 The Interaction Term 


Remark that the interaction in the Hamiltonian is not of the usual form 2 “^) with 

rL/2 

Nl = / dx a*{x)a{x) (2.3) 

J-L/2 

fceIN 

but of the form where 

00 

Nl = y^^ajak- 

k=2 

The reason for choosing Nl in stead of iVi is the breaking of the spatial translation in¬ 
variance in the two terms with A; = 0 and A; = 1. Moreover the latter two terms yield also 
gauge symmetry breaking under the effect of space translations. Using the straightforward 
computation 

cosh(|(T|(a; -f a)) = cosh(|(T|a;)e“l'^'“ -|- sinh(|(T|a) 
for a G IR, we get (for large L): 

n{(^q) ~ aQe“^‘^^“ + (oq — Oi) sinh(|(T|a) 

Ta{ci-o) ~ aoe“l°'l“ -|- (oo — sinh(|(T|a) 

where Tq is the translation automorphism over the distance a G IR. One gets a similar 
expression for Ta{a{). 

In particular, the terms Ta(iVo) and Ta{Ni) are not translation invariant and diverge expo¬ 
nentially for |a| —cx) as 

Taia^ao) ^ 

moreover, these terms break the gauge symmetry. 

Remark also that on the other hand the particle number operator Nl is a good local 
approximation of Nl (2.3) for all translation invariant states, such that the interaction 

\ N‘^ 

term is the appropriate mean held term. 


3 Bounds on the correlation function and condensa¬ 
tion 

Our aim is to hnd the equilibrium states of the system in the grand canonical ensemble. The 
equilibrium state at inverse temperature (3 is characterized by the following correlation 
inequality for all L, see [1,2] 

X]) > MX’X) In (3.1) 

Ul[NX*) 
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for all local observables X, where ujl{-) is the grand canonical equilibrium state at chemical 
potential and inverse temperature /5: 

Tr^B,Bexp{-/3i/£_MF} 


with Xl,b the boson Fock space over £^([—L/2, L/2]). 

Concerning the thermodynamic limit {L —> cx)), we perform this limit keeping the total 
density p constant. Therefore the chemical potential pi is now determined by the particle 
density p and is the solution of the particle density equation: for each given density p we 
have 


P = 


(3.2) 


From the correlation inequality (3.1) follows immediately the inequality 


UL {[X\ [Hl^p - plNl.X]] ) > 0. (3.3) 

In this section we focus on the proof of the condensation for the model Hlmf (2-2). In 
this proof we need bounds which we derive from the correlation inequality (3.1) for some 
specihc observables X’s. Due to the special character of the spectrum (2.1), it is necessary 
to distinguish between products of creation and annihilation operators in the 0- or 1- 
mode, and those in the fc-mode (with k > 2). 

The hrst Lemma is valid for all A; G N. 


Lemma 3.1. If j ^ ki (i = 1,..., m), ki ^ ki/ for i ^ i' and m, Ui G Nq for i = 1,... ,m, 
then 


^(^(eUj)-eUkP)^^ (iV,(iVfc, + 

= UB {{N, + l)iN,,riN,,r ■ ■ ■ (3.4) 


Proof: The proof follows from the correlation inequality (3.1) by taking X successively 
equal to 

and 

□ 

For the chemical potential we hnd the same upperbound as in the case of the free Bose 
gas with attractive boundary conditions. 

Lemma 3.2. For A; = 0,1, one has 


and 


Pl < eL{k) 


p 


lim ul < — 

L^OO 


(3.5) 
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Proof: Take X = al where A; = 0,1 in the inequality (3.3). 

The second inequality (3.5) is obtained by taking the thermodynamic limit L tending to 
inhnity. □ 

Lemma 3.3. For each k = 0,1 and fj,L < we have: 

= ehhL{k)-VL) _ 1 

Proof: The result follows from the inequality (3.1) by taking X successively equal to 
Ofc and al with A; = 0,1. □ 

Lemma 3.4. For each n G N and k >2, we have 

Proof: The result follows from the correlation inequality (3.1) with X = ■ □ 

In order to prove condensation, we need a convenient upperbound for uiiNk) for all A; G N. 
This bound is derived in the following Lemma. 

Lemma 3.5. For each k >2 we have 

^ ^cj) _ 1 

where 

Ck{L) = (3 {^L{k) + ^ - o(e“^'''')^ (3.7) 

Proof: By Lemma 3.4 with n = 0: 

^ - 4 "^^+ 4 ) 

From Lemma 3.2 we know that /i^ < —o(e“'^l°'l). It is also easy to see that uJiiXiNk) > 
0. This leads to 

which gives us immediately the result. □ 

Using the results of the previous Lemma’s, we are now ready to prove the existence of 
condensation in the two lowest energy levels. 
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Theorem 3.1. Let pcond be equal to 


Pcond = lim -ul (iVo + Ni). 

L—5-00 L/ 

Then 

I foo ^ 

Pcond > p - - ^ (3-8) 

where Pcond is the density of the condensate. The condensate density is localized in the 2 
lowest energy levels. 

Proof: From the definition of the particle density p (3.2), we have 

1 1 °° 

—(iVo + Ni) ~ P ~ 

k=2 

By using the estimate of Lemma 3.5, one gets 

1 1 1 

k=2 

with the Cfc(L)’s as in (3.7). 

Taking the thermodynamic limit L —cx) gives us the result (3.8). □ 


Clearly (3.8) shows condensation. Indeed, remark that the integral is convergent for all 
a 7 ^ 0 and that it decreases for (3 increasing. Hence for p large enough or for (3 large 
enough, it follows that the condensate density pcond is strictly positive. 

Finally we derive a result about the type of condensation. We prove that the conden¬ 
sate density is realized in both the two lowest energy modes, with equal weight in the 
thermodynamic limit. 


Theorem 3.2. (i) The condensate is equally distributed on the two lowest energy levels, 

(a) From this, one can compute the asymptotics of the chemical potential p,L for large L: 


P-L = - 0 -^ 


2 


fdpcondL 


+ o{L ^) 


(3.9) 


Proof: 


(i) From Lemma 3.1, with m = 1, j = 1 and k = 0, we get 

ujl{NiNo) + = ool{No) 

Using the spectral properties 

6l( 0) = - 0(e-^''^!) 

6^(1) = -a2 + 0(e-^'"') 
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then 


L L ^ ^ 

Taking the thermodynamic limit L ^ oo, one gets (i). 

(ii) From Lemma 3.3 and since we know that the condensate is equally distributed on 
the two lowest energy levels, we have 

1 . _ 1 . 1 1 _ Pcond 

L “ Le/5hi(0)-/^^) - 1 “ 

Series expansion of this expression with respect to the quantities 6^(0) — hl and 
€ 2 ,( 0 ) = —— 0(e“^l'^l), gives us the asymptotics of p/,. 

□ 
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